Holographic laser Doppler ophthalmoscopy. 
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We report laser Doppler ophthalmoscopic fundus imaging in the rat eye, with near infrared het- 
erodyne holography. Sequential sampling of the beat of the reflected radiation against a frequency- 
shifted optical local oscillator is made onto an array detector. Wide- field maps of fluctuation spectra 
in the 10 Hz to 25 kHz band exhibit angiographic contrasts in the retinal vascular tree without re- 
quirement of exogenous marker. 
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Retinal blood flow plays a pivotal role in several blind- 
ing diseases such as diabetic retinopathy and vascular 
occlusions. While currently available optical instrumen- 
tation is well adapted for imaging retinal vessels, there 
are still technological limitations of currently available 
methods for retinal blood flow measurement that impair 
their clinical applications. Self-beating intensity fluctu- 
ations from a laser spot focused in the eye fundus was 
demonstrated to enable retinal flow assessment [l|-[3| non- 
invasively. It paved the way for the development of spa- 
tial scanning techniques such as confocal spot [2] or line 
jsf scanning Doppler ophthalmoscopy. Yet, laser Doppler 
methods for retinal blood flow mapping are still limited 
in term of velocity resolution and mapping capabilities. 

Optical coherence tomography, commonly used for 
structural retinal layers imaging, was translated to a 
depth-resolved functional Doppler-contrast technique [gI, 
0]- Optical micro- angiography further improved depth- 
resolved Doppler measurements for eye fundus imaging 
(sf. Such techniques can enable quantitative and direc- 
tional flow assessment, but vessel segmentation is re- 
quired [i,[io|. 

Our approach to functional Doppler imaging is based 
on the detection of light fluctuations with an array detec- 
tor. Multiply-scattered light yields Doppler spectra from 
which the directional information of flow is lost because 
of wave vector randomization. Nevertheless, it provides a 
spatially-resolved hemodynamic contrast in low-light il- 
lumination conditions. In this letter, we demonstrate ex- 
perimentally the feasibility of wide-field holographic laser 
Doppler ophthalmoscopy, in vivo. 

The experimental ophthalmoscope realized for this 
study is based on the heterodyne imaging scheme de- 
scribed in [ll|. It consists of a Mach-Zehnder laser in- 
terferometer in off-axis and frequency-shifting configura- 
tion. The detection scheme is sketched in Fig(TJ A laser 
diode provides the main near-infrared radiation at wave- 
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FIG. 1: Optical configuration. The eye fundus is illuminated 
in broad field through the dilated iris. The cornea curvature 
is compensated by a coverslip and ophthalmic aqueous gel. 



length A = 785 nm, polarized linearly. In the object arm, 
a polarizing beam splitter cube is used to illuminate the 
preparation and collect the cross-polarized backscattered 
light component, in order to increase the relative weight 
of multiply scattered Doppler- shifted photons with re- 
spect to photons scattered once 



12|. Three adult rats 



were used for the preparations. Anesthesia was induced 
by intraperitoneal injection of 100 mg/kg ketamine and 
25 mg/kg xylazine (both from Sigma- Aldrich). Topical 
tropicamide (CibaVision) was administered for pupil di- 
lation. Each rat was placed on its side under the illumi- 
nation beam. The head was supported so that the iris 
was perpendicular to the illumination axis. After admin- 
istration of topical oxybuprocaine (CibaVision), a cover- 
slip was applied on a ring surrounding the globe. Methyl- 
cellulose (Goniosol) was applied as contact medium. The 
incident optical field E\ is expanded to form a plane wave. 
Illumination power within 5x5 mm^ at the object plane 
is ~ 1 mW. In the reference arm (local oscillator : LO), 
an attenuator, a half wave plate, and a beam expander 
(not shown) are used to control the beam power, polar- 
ization angle, and to ensure a roughly flat LO illumina- 
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tion of the detector. The optical frequency detuning Au 
between the two optical channels is imposed by acousto- 
optic modulators. The backscattered field E is combined 
with the LO field £^lo with a non-polarizing beam split- 
ter cube. The detuning Auj shifts a component of inter- 
est of the scattered field temporal fluctuation spectrum 
within the actual camera bandwidth (temporal hetero- 
dyning). Moreover, a small angular tilt 6>o,0o of ~ 1° 
ensures off-axis mixing conditions that shift the spatial 
frequency spectrum of the recorded object field (spatial 
heterodyning). The interference pattern I = \E -\- Ei^of 
is measured by a Sony ICX 285AL charge-coupled device 
(CCD) array sensor (gain: 3.8 e/count), from which the 
central 1024 x 1024 pixels region is readout at 10 Hz. The 
detector is set ~ 30 cm away from the object plane. The 
recorded intensity pattern In at time tn in the detector 
plane {z = zq) is : 



In = i^r + i^Lor 



(1) 



where S and £^lo are the complex envelopes of the scat- 
tered and LO field, respectively. Akx = 27rsin(^o)/A, 
Aky = 27rsin(^o)/A are the projections of the difference 
between the LO and signal wave vectors in the transverse 
(x, y) plane. * denotes the complex conjugate. The two 
first terms of the right member of eq. [T] are the self- 
beating (homodyne) contributions of E and £^lo- The 
heterodyne signal of interest lies is the third term. The 
fourth term is the twin-image (ghost) contribution. 

A typical interference pattern In of E beating against 
£^LO is shown in fig. [2fa). The magnitude of the fast 
spatial Fourier transform (FFT) of one recorded frame 
|FFT(/n)P mainly carries the LO contribution in high 
heterodyne gain regime, when the optical power in the 
LO channel is larger than in the object channel [fig. EJb)]. 
Since the LO beam exhibits flat-field fluctuations, its self- 
beating contribution is gathered in the low frequency re- 
gion (1) in reciprocal space (kx^ky). Making the differ- 
ence of two consecutive frames before (or after) applying 
the spatial Fourier transform yields |FFT(/n+i — /n)P; it 
decreases substantially the relative weight of the LO self- 
beating term, showing up the object-against-LO beating 
term and the ghost term [regions (2) and (4), respectively, 
in fig. Efc)]. The object beam self-beating term also ap- 
pears, in region 3. Detection noise is assessed in region 
(5). Because of off-axis geometry, the heterodyne sig- 
nal contribution is shifted-away, by ±{Akx,Aky), from 
self-beating fields contributions. In on- axis geometry 
{Akx = O^Aky = 0), all interferometric terms would 
overlap spatially. The object field can hence be assessed 
with much better sensitivity than if measured in either 
on-axis-only or off-axis-only conditions and used for nu- 
merical reconstruction of the signal S with a discrete 
Fresnel transform S = |FFT((/n+i - /n)e^^^'^'+^'))P, 
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FIG. 2: Recorded interference pattern; Auj/{27t) — 10 Hz (a). 
FFT of the recorded frame (b). The dominating noise (1) is 
gathered in the low spatial frequencies. FFT of the differ- 
ence of two consecutive frames (c). Heterodyne contributions 
(2,4), homodyne contributions (3), noise (5). Image focus (6) 
by Fresnel transform (d). Magnified view of the retina (e). 
White- light endoscopic view of the retina (f). 



where the distance parameter used for free-space back- 
propagation of the optical field is Az = zi—zq = 26.6 cm. 
Once accurate focus is found, the fundus image appears 
as reported in fig. Efd), in region 6, magnified in fig. 
[2][e). These Doppler images are time-averaged (5); an- 
gular brackets ( ) account for averaging over 32 frames. 
The star-shaped vascular tree of the rat eye fundus is 
clearly visible and consistent with the white-light endo- 
scopic image performed afterwards in the same prepara- 
tion, displayed in fig. EJf). 

Detuning Auj slows down selectively the drift rate of 
the set of fringes associated to a given Doppler compo- 
nent, and sets it within the actual camera bandwidth. 
Frequency-selective eye fundus images {S) of a healthy 
rat at four LO detunings (25 Hz, 845 Hz, 4.6 kHz, and 
22 kHz) are reported in Fig. [3l A contrast reversal is 
observed between vessels and surrounding retinal tissue 
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FIG. 3: Doppler fundus images at four frequency detunings 
Ac(;/(27r). Optical power is displayed in logarithmic scale 
(white is for high signal). [AA^] indicates the measurement 
region of the lines reported in fig. [JJa). Scale bar, 1 mm. 
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FIG. 4: (a) (S) against position along [AA'] of fig. [3l^a), av- 
eraged over 15 pixels in the orthogonal direction. Traces cor- 
respond to logarithmically-spaced detuning frequencies from 
10 Hz to 25 kHz. (b) (S) against Au at positions 1, 2, 3, 4 
from the center of the vessel to its periphery. 



(and most likely in the choroid) from low to high fre- 
quency detunings [figs. [3fa) and (c)]. Fluctuation spec- 
tra lines {S) throughout a vessel cross section are re- 
ported in Fig. m as a function of position (a), and as a 
function of Acj, at four locations within the vessel (b); 
typical spatial standard deviation of {S) are shown in 
the inset. It can be estimated that choroidal flow signif- 
icantly contributed to the signal along the retinal vessel 
section [^^']. Within this vessel, the Doppler spectrum 
is clearly broader from those of immediate surrounding 



areas, showing that the retinal flow predominantly con- 
tributed to the signal in this specific zone. The broadest 
Doppler signal, which is still visible at 22 kHz in fig. 
md) is found in the optic nerve head region (circled); it 
is interpreted as a consequence of vessels orientation and 
increased density in this area. 

In conclusion, we have demonstrated the feasibility of 
holographic laser Doppler ophthalmoscopy with near in- 
frared radiation. The illumination power over the whole 
eye fundus can be kept low 1 mW). High detection 
sensitivity, in moderate to high heterodyne gain regime, 
is achieved by spatiotemporal heterodyning, which en- 
ables an efficient rejection of parasitic interferometric 
contributions. Fluctuation spectra discrepancies below 
25 kHz between the superficially vascularized regions and 
the surrounding retinal tissue layers provide an optical 
contrast suitable for angiographic mapping. A high re- 
producibility of the signal acquired over up to several 
hours, over repeated trials, and between different animals 
is observed. Potential applications are essentially the in- 
vestigation of retinal and possibly choroidal vascular dis- 
eases. At the current stage there are yet hmitations due 
to suboptimal lateral, depth, and time resolution, which 
for instance do not allow heart-beat related flow varia- 
tions detection. Technical improvements are expected to 
circumvent these problems. 
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